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Abstract. 

Combining tree decomposition and transfer matrix techniques provides a very 
general algorithm for computing exact partition functions of statistical models defined 
on arbitrary graphs. The algorithm is particularly efficient in the case of planar 
graphs. We illustrate it by computing the Potts model partition functions and 
chromatic polynomials (the number of proper vertex colourings using Q colours) for 
large samples of random planar graphs with up to iV = 100 vertices. In the latter case, 
our algorithm yields a sub-exponential average running time of ^ exp(1.516\/iV'), a 
substantial improvement over the exponential running time ~ exp(0.245A^) provided 
by the hitherto best known algorithm. We study the statistics of chromatic roots of 
random planar graphs in some detail, comparing the findings with results for finite 
pieces of a regular lattice. 

PACS numbers: 02.70.-c, 05.50-f-q, 89.70.Eg 
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1. Introduction 

A typical problem in statistical physics is to compute the partition function Z of some 
model with short-ranged interactions between discrete degrees of freedom defined on 
the vertices of some graph. The partition function is a weighted sum over the states 
of these degrees of freedom. If the graph is a regular lattice (i.e., consists of a large 
number M of identical layers) one usually rewrites Z in terms of matrix elements of 
T*^, where T is a transfer matrix (or imaginary-time evolution operator in an equivalent 
path integral formulation) that corresponds to the addition of a single layer to the 
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lattice. This suggests solving the problem by diagonalising T, so that, in some sense, 
one has substituted algebraic complexity for combinatorial complexity. When the lattice 
is planar, the exact diagonalisation of T in the limit of an infinitely large system (or 
thermodynamic limit) can in many cases be achieved by using the powerful tools of 
integrability [Ij. Alternatively, there exists many efficient means of diagonalising T 
numerically for rather large systems. 

In many applications one is however interested in models defined on graphs that are 
not regular, but incorporate some element of randomness. One must then distinguish 
whether this randomness is of the annealed or the quenched type. For the annealed 
disorder, Z is a double sum over the graphs and the vertex degrees of freedom. In the 
case of planar graphs, this double sum can in many cases be evaluated analytically by 
random matrix techniques [2]. But arguably, the physically most relevant scenario is 
that of quenched disorder, where one would typically need to average the free energy 
logZ, not just Z, over the disorder realisations. 

Clearly, a maximal way to deal with quenched randomness would be to evaluate Z 
independently for each graph in the sample. This is generally not possible analytically, 
even for planar graphs. The purpose of this paper is to exhibit an efficient algorithm that 
evaluates Z exactly for an arbitrarily given graph. This algorithm applies to a rather 
large class of models of statistical physics, as outlined above, and it does not require 
the graph to be planar. Rather than delving into this generality — we shall briefly come 
back to the issue in the discussion — we have here chosen to focus on one signiflcant 
speciflc application, namely the evaluation of the Potts model [3] partition function 
Zg{Q,v) (which, after a change of variables, equals the Tutte polynomial [1] known in 
graph theory) for a given planar graph G. We are particularly interested in the special 
case V = —1, where Zg{Q,v) = Xg{Q) equals the number of Q-colourings (chromatic 
polynomial) of G. We now outline the physical motivation for studying this problem, 
before ending the introduction by giving a brief account on the working principle of our 
algorithm and its performance (computational complexity). 

The Q-colouring problem consists in assigning to each of the vertices of a graph 
G any one of Q different colours, in such a way that adjacent vertices carry different 
colours. Each such assignment is known as a proper vertex colouring. The Q-colouring 
problem arises within physics in studies of frustrated antiferromagnetic spin models, and 
in spin glass theory in particular [5]. The number of possible colourings (possibly zero) 
can be shown [Gj [7] to be a polynomial in Q, known as the chromatic polynomial Xg{Q) 
of the graph. In particular it makes sense — and is useful — to generalise the original 
counting problem and to study the properties of Xg{Q), considered as a polynomial of 
a formal variable Q. 

The history of the Q-colouring problem is long and interesting, and we refer the 
reader to [8] for an extensive list of references. The case where G is a planar graph has 
attracted particular interest. A well-known result is then the four-colouring theorem 
[9] which states that Xg(4) > for any planar G. Other results exploit the extension 
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of Xg{Q) to non- integer values of Q. One interesting question is whether there exists, 
in the planar case, some Qc so that Xg{Q) > for all Q G [Qc,oo). This statement 
has been established as a theorem [10] for Qc = 5, and the so-called Birkhoff-Lewis 
conjecture [10] that this extends to Qc = 4 is widely beheved to be true. 

The case where G is a regular planar lattice has also been intensively studied. 
In particular, the location and properties of chromatic roots, Xg{Q) = 0, in the 
complex Q-plane has been studied for a variety of lattices and boundary conditions 
(see [m [HI 1121 1131 El US] and references therein). Some of the mechanisms responsible 
for the generation of real chromatic roots in the region Q G [0, 4], close to or exactly at 
the so-called Beraha numbers Bk = {2 cos{tt / k))'^ with k > 2 integer, have even been 
understood analytically [161 112] • One should also mention that there exists a family of 
planar graphs — planar triangulations actually — with real chromatic roots converging to 
Q = 4 from below p^, meaning that the value of Qc in the Birkhoff-Lewis conjecture 
cannot be lowered further. 

One issue which has been left largely unanswered by these studies is the distribution 
of the location of (real and complex) chromatic roots for ensembles of randomly chosen 
planar graphs. This situation appears to be the most interesting from the point of view 
of the physics of disordered frustrated systems. From the mathematical side, it has been 
proven [19] that chromatic roots are dense in the complex plane (except maybe in the 
disc IQ — 1| < 1) for a special class of planar graphs (generalised theta graphs). However, 
this might not have much to do with the physical question of the chromatic roots of a 
typical planar graph. One would also want to know whether typical roots accumulate at 
the Beraha numbers, and which graph characteristics (connectivity, average coordination 
number,. . . ) might be responsible for the location of roots. 

In order to elucidate the amazingly intricate behaviour of chromatic roots in the 
limit of large graphs, it is clearly desirable to have efficient means of computing Xg{Q)- 
Let us recall that in theoretical computer science, an important outstanding question 
is whether the class P of decision problems that can be answered in polynomial time 
coincides with the class NP of problems for which a proposed answer can be verified 
in polynomial time. NP-complete problems are those to which any problem in NP can 
be reduced in polynomial time. At present, no polynomial-time algorithm has been 
found for any of the thousands of known NP-complete problems and it is hence widely 
believed that P ^ NP. Likewise, one can define a counting analogue of NP, denoted 
by #P, as the class of enumeration problems in which the objects being counted are 
the possible answers accepted by an NP machine. If an NP problem is often of the 
form "Are there any solutions that satisfy certain constraints?" , the corresponding ^P 
problems ask "how many" rather than "are there any". The class T^P-complete of the 
hardest problems in t^P is defined likewise. Clearly, #P-complete problems are at least 
as hard as NP-complete problems. It is known [20] that the counting of proper vertex 
colourings is a #P-complete problem for Q = 3, and the same is then true a fortiori for 
the computation of Xg{Q), let alone Zg{Q,v) [21j, for generic values of Q and v. 
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In practice this means that any algorithm computing Xg{Q) can be expected to have 
a running time that increases exponentially with the number of vertices A^. However, 
lowering the coefficient of the exponent can still make a huge difference for studying the 
issues outlined above. One central goal of this paper is to make a substantial step in 
that direction, by lowering the average running time ~ exp(0.245iV) of the previously 
best known algorithm J22j to ~ exp(1.516A/iV), as illustrated in Fig. [1] below. The 
improvement from exponential to sub-exponential asymptotics is not only important 
from a theoretical point of view. To get a rough idea of what this means in practice, in 
ten seconds the algorithm [22] will compute Xg{Q) for & typical planar graph of = 40 
vertices, whereas our algorithm can deal with = 100 in the same time. 

Algorithmic progress on T^P-complete problems related to graph theory and 
network design has been made by several, usually widely separated, communities. 

On one hand, statistical physicists have shown that the relevant partition functions 
can be constructed in analogy with the path integral formulation of quantum mechanics. 
To this end, the configuration of a partially elaborated graph are encoded as suitable 
quantum states, and the constant-time surface is swept over the graph by means of 
a time evolution operator known as the transfer matrix. Although rarely stated, this 
approach is valid not only for regular lattices but also for arbitrary graphs. 

On the other hand, graph theorists have used that graphs can be divided into 
"weakly interacting" subgraphs through a so-called tree decomposition [23], El], and 
solutions obtained for the subgraphs can be recursively combined into a complete 
solution. 

The principle underlying the algorithmic progress to be described in this paper 
is that tree decomposition and transfer matrix methods can be combined in a very 
natural way. The main idea is that the tree decomposition is compatible with a recursive 
generalisation of the time evolution concept. Borrowing ideas from quantum field theory, 
the combination of partial solutions is obtained by the fusion of suitable state spaces. 
The resulting algorithm works on any graph, and can readily be adapted to many other 
problems of statistical physics, by suitable modifications of the state spaces and the 
fusion procedure. 

In particular, we will apply this technique to the problem of computing the 
chromatic polynomial on planar graphs, obtaining exact solutions, for graphs with 
~ 100 vertices, in only a few seconds. 

The layout of the paper is as follows. In section [2] we briefly recall the relation 
between the Potts model partition function Zg{Q,v) and the chromatic polynomial 
Xg{Q)- III section [3] we present our algorithm and discuss its performance. The results 
for chromatic roots of planar graphs are given and discussed in section |H Finally, we 
give our conclusions in section |5l 
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2. Potts model and vertex colourings 

We first recall the relation between the vertex colouring problem and the Potts model. 
Consider a graph G = (V, E) with vertices V and edges E, and let at = 1,2, ... ,Q be 
the colour of vertex i & V. Then 

= E n ^'"'^'"'''''^ (1) 

a {ij)eE 

is the partition function of the Potts model on G. The Kronecker delta 6{ai, aj) = 1 if 
(Ti = aj, and otherwise. Inserting the obvious identity e^^^"^^'"^^ = 1 + vS{ai, aj), with 
f = e^ — 1, and expanding out the product we obtain the Fortuin-Kasteleyn expansion 
[25] 



^G(g,^) = E^'^'^'^^^ (2) 

ACE 

where k{A) is the number of connected components in the subgraph G' = {V,A). 
Obviously, in the antiferromagnetic limit K — )■ — oo (or equivalently v — —1) the only 
surviving configurations are proper Q-colourings of the graph G. Indeed, the special 
case Xg{Q) = Zg{Q, —1) is a polynomial in Q, known as the chromatic polynomial, and 
equals the number of vertex colourings. 

3. Algorithm 

3.1. Transfer matrix 

The computation of the partition function Zg{Q,v) using the expansion ([21) has been 
described in [26l [8] for the case where G is a finite piece of a regular lattice (notably, 
but not exclusively, a planar one). 

We first describe how this traditional transfer matrix method extends to the 
case where G is an arbitrary graph — i.e., not necessarily part of a regular lattice, 
nor necessarily planar. In short, the combined action of linear operators builds 
a superposition of all configurations appearing in the partition funciton with their 
correct statistical (Boltzmann) weight entering as coefficients. To better illustrate this 
procedure, consider the following example graph G 

6 

/ ^7 

1 — 3 ^ / 

(3) 

2 — 4 ^ \ 

\ .9 
8 

We first have to define the order {vt} in which vertices will be processed. This order is 
the basis for the construction of a "time slicing" of the graph. With each time step is 
associated a bag (a vertex subset) of active vertices. A vertex becomes active as soon as 
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one of its neighbours is processed and it stays active until it is processed itself. Taking 
the vertices in lexicographic order we obtain the following decomposition: 



(4) 



where we wrote in bold face the vertex being processed at each time step. 

Each bag has its own set of basis states consisting of the partition^ of the currently 
active vertices. For instance, in the first time step, the basis states are the five partitions 
of the three-element set {1, 2, 3}: 

I o o o \ I o o o \ I o o o \ I ©""cT^o \ I o o o \ 
Il23'' Il23'' Il23'' Il23'' Il23'' 

These partitions describe how the active vertices are interconnected through A (1 Et, 
where Et ^ E is the subset of edges having been processed at time t. A state is a linear 
superposition of basis states. 

Processing a vertex consists in processing edges connecting it to unprocessed 
vertices and then deleting it. Since each edge e E E may or may not be present in 
A we process an edge by acting on the state with an operator of the form 1 + f J^- 
where 1 is the identity operator and Jj-,- a join operator. A join operator acts on a basis 
state by amalgamating the blocks containing vertices i and j. 

Vertex deletion is defined in terms of a deletion operator Dj that removes i from the 
partiton and applies a factor Q (resp. 1) if i was (resp. was not) a singleton. 

D,|o o ...) = g|o ...), (6) 

D.lfo ...) = |o ...). (7) 

For example, in (j4]), processing the first bag means computing the following composition 
Di(l + t;Ji2)(l + t;Ji3) I ° ° °), 



which gives 



{Q + 2v)\o o )+v^\o o) 



2 3' 123 

concluding the first time step. 

When a new active vertex is encountered it is inserted, as a singleton, in each 
partition composing the current state. After processing the last bag, the complete 
partition function ([2]) is obtained as the coefficient of the empty partition resulting from 
the deletion of the last active vertex. 

At each step, the time and memory requirements are determined by the bag size n 
which is the number of vertices simultaneously active. If the graph is planar, the number 

§ Recall that a partition of a finite set is a collection of nonempty, pairwise disjoint subsets of S 
whose union is S. The subsets are called blocks. 
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of partitions to be considered is at most the Catalan number C„, i.e., the number of 
non-crossing partitions. The corresponding generating function is 



C(.) = ^C„- = i^l^|Hl. (8) 

n=0 

If the graph is not planar, the number of partitions is at most the Bell number Bn, with 
generating function 



i?(z) = Vi?„-=exp(e^-l). (9) 

n=0 

We have Cn = 4"n^^/^7r^-'^/^[l + 0(l/n)], whereas i?„ grows super-exponentiallylj] 



3.2. Tree decomposition 

It turns out that the decomposition (jlj) of G is a special case of a more general 
construction. By definition, a tree decomposition [23l [24] of a graph G = (V, E) is a 
collection of bags, organised as a tree (a connected graph with no cycles), and satisfying 
the following requirements: 

(i) For each i & V, there exists a bag containing i; 

(ii) For each (ij) G E, there exists a bag containing both i and j; 

(iii) For any i & V, the set of bags containing i is connected in the tree. 

The previous decomposition @ is just a special case of a tree decomposition (a path 
decomposition) . As an example of the general construction, applied to ([3]), consider 

where the arrows form the unique path that connects each bag to the central one (the 
root of the tree). 

The advantage of working with tree instead of path decompositions relies on the 
fact that in the former case a decomposition with smaller bags can be obtained (the 
latter being just a special case). Therefore, the number of states one has to keep track 
of is exponentially smaller, and the gain is significant. 

The transfer matrix approach can be adapted naturally to this new general setting: 
Properties (i)-(ii) guarantee that each edge and vertex are processed within a definite 
bag. Property (iii) implies that each vertex has a definite life time in the recursion, its 
insertion and deletion being separated by the processing of all edges incident on it. 

I Note that the restriction to non-crossing partitions holds true even if the ordering of the vertices, 
such as in ([4]), does not respect the planar embedding ([3|) of the whole graph (which will in general 
be unknown). Because we choose to store partitions independently of their crossing property, it has 
not been relevant to be able to explicitly keep track of the planar embedding. All that matters is that 
the number of partitions does not exceed C„, and this is true since such a non-crossing representation 
could be obtained "by hand" by inspecting the planar embedding at each stage of the transfer process. 
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In this new version the algorithm starts from the root of the tree, which can be 
chosen arbitrarily, and runs through the tree recursively. Alternatively, this can be 
thought of as starting from the leaves and building up the tree inductively. 

Consider first the case of a parent bag P with only one daughter bag D. Going 
up from D to P in the recursive step involves deleting vertices D\P, inserting vertices 
P \ D and finally processing edges in P. A tree decomposition does not specify when 
an edge e = (ij) must be processed. A simple recipe would be to process e as soon as 
one encounters a bag containing both i and j. However we note that this freedom of 
choice can be exploited to optimise the algorithm (see below). 

3.3. Fusion 

We now discuss the case when a parent bag P has several daughters Di with £ = 
1,2, ... ,d. In this case, vertex deletions and insertions are followed by a special fusion 
procedure. Suppose first d = 2, and let Vi be a partition of fl P with weight wi, and 
V2 a partition of H P with weight W2. The fused state is then 



and it occurs with weight W1W2. Here V denotes the join operation in the partition 
latticefl For definiteness, let us explain how this can be computed in practice. First 
choose some set Ei so that 



where Si is the all-singleton partition of Di fl P. Note that it is not necessary to choose 
El as a subset of E, although this can always be done. Since JgJe' = Je'Je the order in 
the above product is irrelevant. Similarly choose E2 for 7^2- The fused state can then 
be constructed explicitly as 



where Su is the all-singleton partition of {Di U D2) fl P. For d > 2 daughters, the 
complete fusion can be accomplished by fusing Di with D2, then fusing the result with 
-D3, and so on. 

% We recall that a partial order of the partitions follows from defining a partition Va to be a refinement 
of Vb, and we write Va 'Pb, provided that each block in Va is a subset of some block in Vb- With 
this partial order, the set of partitions form a lattice, in the sense that any two partitions Vi and V2 
possess a largest lower bound called meet and denoted Vi A V2 and a smallest upper bound called join 
and denoted Vi V 7-'2 ■ The meet Vi A V2 has as blocks all nonempty intersections of a block from Vi 
with a block from 7^2- The blocks of the join Vi\/V2 are the smallest subsets which are exactly a union 
of blocks from both Vi and 7^2- In the partition lattice, the bottom (or finest) element is the unique 
partition in which each block has size 1 (the all-singleton partition) , and the top (or coarsest) element 
is the unique partition in which there is a single block (the all-connected partition). 



Vi®V2 = Viy V2 




(11) 




12 , 
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Let us illustrate the fusion procedure for G with the tree decomposition (ITU]) . After 
processing the two left-most bags and deleting vertex 2, the propagating state is 

^^ll °3 2)+^2| fo) = (Wl+W2J34)| O O), (12) 

where Ui = + 3v{Q + v) and U2 = Q'^v + ?>Qv^ + Av^ + f ^. By symmetry, the same 
result is obtained for the two top right bags (replacing 4 by 5) and for the two bottom 
right bags (replacing 3 by 5)|l| The fused state arriving in the central bag is then 

= ^1 I g I %) + (3t^it^2 + ^^2) I %Y\ ) 

from which the result Zg{Q,v) follows upon deleting vertices 3,4,5. 
3.4- Pruning 

Problem specific features can be exploited to reduce further the number of basis states 
to be considered. As an example of this, note that in the colouring case (f = — 1) the 
operator Ojj = 1 + v^ij associated with an edge (ij) is a projector, 0^^- = Oj^, and it 
annihilates the subspace of partitions where i and j are connected. It follows that one 
can discard basis states in which two vertices are connected, as soon as one discovers 
that an edge between them is going to be processed later within the same bag or within 
the parent bag. Especially before fusions this simple trick reduces substantially the 
number of basis states and thus leads to a big speed up. 

3.5. Performance 

For a planar graph, the state of a bag of size n is spanned by C„ basis states. (For a 
non-planar graph, replace C„ by The memory needed by the algorithm is therefore 
proportional to C„^^^, where rimax is the size of the largest bag. The time needed to 
process one edge in a bag of size n is proportional to C„. 

However, most of the time is spent fusing states. For a parent P with d daughters 
Di, the number of basis state pairs to be fused is 

d 

^C\v,,^nP\C\DtnP\, (14) 
e=i 

where we have set Vk = U'l^De. Each of these elementary fusions (i.e., the joins Vi VP2) 
can be done in time lineaiu in the number of participating vertices. Note that we can 
choose the order of successive fusions so as to minimise the quantity (fT^ . 

+ Needless to say, the fusion procedure will work also in cases where the graph possesses no such 
symmetries. In our example, we have chosen C? as a rather symmetric graph in order to keep the actual 
computations as simple as possible. 

* We represent partitions as lists of numbers linking each participating vertex to the block it belongs 
to. The set Ei in (fTTj) can be obtained in a single scan of this list where one keeps track of the last 



Tree-decomposed transfer matrix 



10 



It is therefore essential for the algorithm that one knows how to obtain a good 
tree decomposition. The minimum of n^ax ~ 1 over all tree decompositions is known as 
the tree width k, but obtaining this is another NP-hard problem. However, the simple 
algorithm GreedyFillln [27] gives an upper bound ko on k and a tree decomposition 
of width ko in time linear in the number of vertices N. For uniformly generated planar 
graphs we find that for N = 40 — a value enabling comparison with algorithms that 
determine k exactly — that k^ = k nearly always {ko = k + 1 with probability ~ 10^"^). 
When ko is small enough that all the partitions can fit into the computer's memory, the 
algorithm has proven to be very fast with execution time in the order of seconds. 

To summarise, suppose that we wish to compute Zg{Q,v) or Xg{Q) for a planar 
graph of vertices and M edges, and that we are given a tree decomposition of B 
bags, with rzmax being the size of the largest bag. The number of (binary) fusions to 
be performed is F = ^j(c?i — 1), where di is the number of daughters of bag i. For 
simplicity we assume a computational model where the operations on the weights can be 
done in unit timej^ For the version of the algorithm without pruning, we have therefore 
shown that the worst-case running time is asymptotically proportional to 

{N + M + B)Cn^^^n^.. + i^C^Lax^max • (15) 

Note that this is linear in N and M for a fixed nmax (cf. [28]). 

We choose to test our algorithm against the one presented by Haggard et al. in 
[22]. We first generated a uniform sample of 100 planar graphs for each size = 
between 20 and 100 using Fusy's algorithm [31]. We then ran four different algorithms 
over this sample: the algorithm of [22], our first path-based transfer matrix algorithm, 
the new tree-based version algorithm and a tree-based version using the above pruning 
optimisation. Path decompositions were obtained with a variant of the GreedyFillln 
algorithm in which the resulting tree decomposition was forced not to branch. Average 
running times are presented in Fig. [H in the form of effective exponential fits. 

While these fits clearly confirm the efficiency of the tree decomposed approach, 
they are actually misleading and hide an important fact. Namely, the tree width 
of a planar graph is bounded by a\/iV, as follows from the famous planar separator 
theorem [29]. The currently best upper bound on the constant is a < 3.182 [30j. These 
facts — combined with (|T5l) and ([8]), and the heuristic observation of the near-ideality 
of the tree decompositions obtained by the GreedyFillln heuristics — immediately 
imply an upper bound on the worst-case running time on the tree-based algorithms of 
]^g3.i82v^ _ exp(8.822\/iV). It follows that the mean running time of these algorithms 
must be of the form exp(/3-\/iV), with the constants (5 satisfying /3troe+pruning < Aree- The 
fits for /3, using the same data as in Fig. [H are shown in Fig. [2j 

seen vertex for each block. In this representation, applying Jg is constant time, apart from a possible 
linear-time procedure to bring the resulting list in a canonical form. This last procedure can be safely 
postponed after all the join operations, making the whole elementary fusion linear in the number of 
participating vertices. [Note that in [28] the join operation is claimed to be quadratic in n] 
tt The weights are numbers for an evaluation of Zq{Q,v) or of xg{Q)^ polynomials with integer 
coefficients for xg(Q), and polynomials in two variables with integer coefficients for Zg{Q,v). 
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Figure 1. Average running time in seconds on a random planar graph of size N 
vertices. Each point is averaged over 100 graphs. 



3.6. Technical aspects 

For completeness we give a few implementational details. Partitions relevant to the 
propagating state are kept in a hash table for fast (amortised constant time) access. 
The corresponding weights are polynomials in Q whose coefficients rapidly exceed the 
machine integer size M = when running on big graphs. To solve this issue without 
requiring additional memory we used modular arithmetics: the algorithm was run many 
times with coefficients computed modulo primes p < M, and the original coefficients 
were reconstructed by the Chinese remainder theorem. 

4. Results 

As a simple application of our algorithm we obtained the distribution of the chromatic 
roots {Q I Xg{Q) = 0} for ensembles of planar graphs of up to = 100 vertices. This 
problem is interesting both in statistical mechanics and in graph theory [U [161 [IZl E]- 
As known from Lee- Yang theory, the roots signal phase transitions. 

Before turning to the actual results, it is useful to give some results for regular 
planar graphs of a few hundred vertices. To this end, we have used a particular 
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Figure 2. Average running time in seconds on a random planar graph of size N 
vertices, now fitted to the form exp(/3\/^)- The data points are identical to those of 

Fig.m 

adaptation of the traditional (path-decomposed) transfer matrix algorithm that takes 
special advantage of the regular lattice structure ^32j to compute Xg{Q) for LxL sections 
of the square and triangular lattice (the latter being considered as a square lattice with 
added diagonals). The boundary conditions are free-free, in the terminology of [8]. The 
chromatic polynomials for L = 10, 12 have been checked against [12], while those for 
L = 14, 16, 18 are new. 

The resulting locations of the chromatic roots in the complex Q-plane are shown 
in Fig. |3] for the square lattice, and in Fig. H] for the triangular lattice. As in earlier 
work [SI [121 US] we notice that most of the roots fall on connected curves which, very 
roughly speaking, tend to form a egg-shaped figure with a pair of prongs on the right 
side. Without going into details (which are numerous, see [331 HB ESI [13 El [111 [131 [M] 
for more exact statements) in the limit L — )• oo the egg-shaped part is expected to 
enclose a segment Q G [0, Qq] of the real axis, with Qq = 3 for the square lattice, and 
Qo ~ 3.81967 for the triangular lattice. In the interior of the egg one observes additional 
isolated real zeros right at, or extremely close to, the so-called Beraha numbers 

Bk = i2cos{n/k)f (fc = 2,3,4,...) (16) 




U I I I I I I I I I I I I u 

0.5 1 1.5 2 2.5 3 
Re(Q) 

Figure 3. Chromatic roots for Lx L sections of the square lattice with free boundary 
conditions, and L = 10, 12, 14, 16, 18. Beraha numbers Bk with fc = 2, 3, 4, 5 are shown 
as small vertical line segments. 

The reason why this is so is hnked to particularities of the representation theory of the 
quantum algebra that underlies the two-dimensional Potts model [161 HI]- 

The existence of a finite-size equivalent of Qo is clearly visible from Figs. [3H11 One 
could propose defining Qq{L) as the largest real root, but one should be careful since in 
some cases the egg-shaped curve does not possess a root right on the real axis. Barring 
these (and other) difficulties, one could speculate that if the result for an arbitrary 
planar graph was qualitatively similar to that of these regular lattices, the probability 
distribution of complex roots would look like some broadened-out egg. In addition, the 
distribution of real roots would be a superposition of sharp peaks centered at the Beraha 
numbers, and a broad background distribution corresponding to the Qq{L). 

This is indeed more-or-less what we observe. For a sample of 2500 uniformly drawn 
[3T] planar graphs of = 100 vertices we show the distributions of complex chromatic 
roots in Fig. |5l and the distribution of the subset of real roots in Fig. El Regarding 
the complex roots, we note that although chromatic roots of planar graphs have been 
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Figure 4. Chromatic roots for L x i sections of ttie triangular lattice with free 
boundary conditions, and L = 10, 12, 14, 16, 18. Beraha numbers B}~ with fc = 
2, 3, 4, 5, 6, 7, 8 are shown as small vertical line segments. 

shown to be dense in the complex plane [19] (except maybe in the disc IQ — 1| < 1), 
typical roots are clearly quite close to the origin. 

As to the real chromatic roots, the absense of roots on the negative real axis, and 
the intervals (0, 1) and (1,32/27] follow from a theorem [35] (see also [36])- The roots 
found here respect this theorem as well as the Birkhoff-Lewis conjecture [TU]. Apart 
from that. Fig. |6] shows as expected a superposition of a broad background distribution 
and sharp peaks centered ai Q = with k = 2,3,4,5,6 (we have ~ 2.61803). 
It seems likely that for N ^ oo this background will extend to the interval (32/27,4) 
and peaks will occur at all Beraha numbers. We also expect that the maximum of the 
background distribution may be shifted further to the right by requiring the graphs to 
be two- or three-connected, and so presumably the peaks at the first few Bk should 
stand out more clearly. We plan to investigate this and further issues in more detail 
elsewhere [37] . 



Tree-decomposed transfer matrix 



15 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



X 

Figure 5. Distribution for complex chromatic roots Q = x + iy obtained from a 
sample of 2500 planar graphs of size 100. The normalisation is such that p{x, y) is the 
expected number of roots in [a; — 0.05, a: + 0.05] x [y — 0.05, y + 0.05] counted with their 
multiplicity. 

5. Conclusion 

In this paper we have shown how to obtain an efficient algorithm for solving instances of 
graph-related #P-complete problems. The case of the Potts partition function Zg{Q, v), 
and in particular its specialisation the chromatic polynomial Xg{Q) = Zg{Q, —1), were 
studied in detail, and some results on the distribution of real and complex chromatic 
roots for uniformly drawn random planar graphs were given. 

Our algorithmic approach has some overlap with an algorithm described — but 
apparently not implemented — by Noble [28]. This author also combines partial 
evaluations of Zg{Q,v) — treated by him in the Tutte polynomial formulation [1| — in 
the basis of partitions with the notion of tree decomposition. There are however a 
number of important differences. First, our use of transfer matrices — and specifically of 
its factorisation within each bag, i.e., our choice to process a single edge at a time within 
each bag — leads to a better time complexity. In particular, treating a daughter bag 
with no sisters has worst-case running time ni,^^ + C„ 2"™'''''^"™'"'^-'^)/^ in Noble's 

^ J^max max ' ''-max 
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Figure 6. Probability distribution for the subset of chromatic roots in Fig. [5] that lie 
on the real axis. Normalised to unit total area. 

approach versus C„^^^nmax in ours. Second, we treat binary fusions more efficiently, 
exploiting the fact that a given pair Vi, V2 of input partitions gives rise to a unique 
output partition V\\JV2- this reduces the time complexity for binary fusions from worst- 
case time C'^^^^n^j^^ in [28] to C^^^^rimax in our approach. Third, our algorithm appears 
easier to describe and implement — as we have indeed done — and it avoids having to deal 
with the Q — limit (x = 1 in the notations of [28]) as a particular case. 

It is thus clear that what matters the most is to obtain a good tree decomposition, 
and in particular the running time is essentially determined by rimax- In Fig. [7] we show 
the distribution of rimax obtained from applying the algorithm GreedyFillln to random 
planar graphs as a function of A^. 

It is plainly visible from Fig. [7] that both the mean and worst-case n^^^ exhibit a 
slower than linear growth with A^. Assuming that GreedyFillln yields always an n^ax 
which is of the order of the true tree width, we would have rimax < cst-\/iV, and the data 
of Fig. [7] appear to be compatible with such a scaling. We defer the further investigation 
of the asymptotic behaviour of GreedyFillln to future work [37] . 

We should also mention that the use of tree decomposition in the context of decision 
problems has previously been advocated by a number of authors (see [38] and references 
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Figure 7. Probability distribution p{N, rimax) of obtaining a tree decomposition of 
maximum bag size rimax (i-e., width rimax — 1) when applying the GreedyFillln 
algorithm to random planar graphs, as a function of the graph size N. The data 
is computed over a sample of 100 planar graphs for each size between 20 and 500, 
binned in blocks of 10 on the x axis and normalised to unit total area. 

therein). In these works, the solution is found by doing dynamic programming on the 
tree-decomposed graph. 

Let us conclude by commenting on the generality of our algorithm. Going through 
the list of known NP-complete problems related to graphs and network design, one 
easily realises that most of them can be promoted to counting problems (which cannot 
be easier), and a counting polynomial analogous to Xg{Q) can be defined. Adapting our 
algorithm to those cases usually requires just a problem-specific definition of the states, 
of the single-edge transfer process, and of the fusion procedure, whereas the bulk of the 
method can be taken over without changes. 

It is worth underlining that our ability to solve efficiently an enumeration problem 
can be exploited to obtain a specific instance that solves the corresponding decision 
problem with only an additional linear time factor. 

In particular, we have adapted (but not yet implemented) our algorithm to counting 
versions of the following problems: Hamiltonian walks and cycles, longest path, vertex 
cover, dominating set, feedback vertex set, and minimum maximal independent set. We 
plan to report on these problems elsewhere [37]. 
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